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SUMMARY 
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Numerical solutions of the complete, time -averaged conservation equations using 
several eddy-viscosity models for the Re 3 molds shear stress to close the equations are 
compared with experimental measurements in a compressible, turbulent separated flow. 
An efficient time -splitting, explicit difference scheme was used to solve the two- 
dimensional conservation equations. The experiment used for comparison was a turbu- 
lent boimdary layer that was separated by an incident shock wave in a Mach 2.93 flow 
with a unit Reynolds number of 5.7 x lO'^/m. Comparisons of predicted and experimental 
values of surface pressure, shear stress along the wall, and velocity profiles are shown. 
One of the tested eddy-viscosity models which allows the shear stress to be out of equi- 
librium with the mean flow produces substantially better agreement with the experimen- 
tal measurements than the simpler models. A tool is thereby provided for inferring 
additional information about the flow, such as static pressures in the stream, which 
might not be directly obtainable from experiments. 


INTRODUCTION 


Development of turbulence models capable of adequately predicting separated com- 
pressible flows requires extensive comparisons of calculations with a variety of experi- 
ments. At present, the proper role of the calculations is perhaps to aid in the design, 
interpretation, and documentation of the experiments. For those purposes, it is expedi- 
ent to seek simple empirical relationships that can be used in the calculations to corre- 
late the experiments. Only after a reliable body of experimental data has been accumu- 
lated will definitive verification of the advanced turbulence models be possible. 

The adequacy of the numerical method for solving the conservation equations has 
been established by comparison of calculations with experimental measurements in 
laminar separated flows (refs. 1 to 3). Methods and procedures for incorporating turbu- 
lence models efficiently in the numerical methods are still in a stage of development. 
Thus, while experience is being gained with these methods and procedures, it is expedi- 
ent to utilize simple empirical relationships in most of the calculations rather than add 
the complication of advanced turbulence models. 
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Several solutions to the conservation equations for compressible turbulent sepa- 
rated flows have been obtained to date (refs. 4 to 8). The present investigation was 
motivated and guided by interaction with experiments conducted by the authors of ref- 
erences 9 to 11. The class of experiments considered in this study (i.e., non -hypersonic, 
near -adiabatic wall condition) was chosen for examination to minimize uncertainty with 
the body of data to be used for comparison with the calculations. In addition, the effects 
of pressure fluctuations which might influence the turbulence modeling are minimized in 
such flows. An attempt has been made in the present study to'fihd modifications of a ^ 
simple algebraic turbulence model that will adequately represent the turbulent shear 
stress throughout a shock -separated turbulent boundary layer. 

SYMBOLS ' ^ ^ 

Cf skin-friction coefficient p<^Uoo^^ ' ■ ' 

Cp pressure coefficient 

IL mixing length 

M Mach number 

p pressure 

Pr Prandtl number (0.72) 

Pr-p turbulent Prandtl number (0.9) 

Re , Reynolds number 

T temperature 

u streamwise velocity 

Umax maximum velocity in boundary layer 

Uj friction velocity 

V velocity in y-direction 


402 



X 


y , . 


streamwise coordinate 

. , J 

. inviscid impingement location ^ 
cross -stream coordinate , . 

.. . ..r (boundary -layer thickness , 


e eddy viscosity 

^Inner eddy viscosity near wall 

^outer eddy viscosity in outer part of boundary layer 

M absolute viscosity 

p density 

Txy shear stress fluxes in Navier -Stokes equations 

T total shear stress, -r^y 

Subscripts: 

e edge of boundary layer 

o upstream conditions 

t total conditions 

w wall conditions 

free -stream conditions 
Svqperscript: 

+ indicates dimensionless boundary-layer quantity 



PRELIMINARY CONSIDERATIONS 


Description of the Flow 

The flow field considered in this study is shown schematically in figure 1, (See 
ref. 10 for complete detail.) The free-stream Mach number is 2.93 and the unit Re 3 molds 
number is 5.7 x lO^/m. A shock wave, generated by a plate set at 13° to the free stream, 
impinges on the turbulent boundary layer at the upper nozzle wall. This incident and 
reflected shock system produces a large pressure rise (Pfinal/Pinitial ~ causes the 

boimdary layer to separate. The pressure rise feeds upstream in the subsonic part of , 
the boundary layer and forms a separation bubble containing reversed flow; hence, sepa- 
ration occurs well upstream of the inviscid shock impingement point. An induced shock 
wave results from the upstream separation and a reflected shock wave also forms. 
Although care was taken in the experiment to minimize three-dimensional effects, the 
flow may still contain these effects as noted in reference 10. This is important to con- 
sider since the present numerical method applies strictly to two-dimensional flows. 


The Numerical Method 

The time-averaged Navier -Stokes equations to be solved are listed in reference 5. 
The numerical method and special numerical procedures used in the calculations pre- 
sented herein are described in detail in references 2 and 5. The grid used in the present 
study had 40 nodes in the streamwise direction with Ax = 5.1 x 10~3 m and 32 nodes in 
the cross -stream direction with four different Ay ranging from 2.5 x 10"® m near the 
wall to 4.8 X 10*^ m beyond the boundary -layer edge. - 

A simplified turbulence model used previously (refs. 5 to 7) to gain experience in 
the numerical solutions while conserving computer time expresses the shear stress as 



( 1 ) 


The eddy viscosity e is modeled in two regions within the boundary layer, the inner 
region comprising essentially the Ic^arithmic region near the waU and the outer region 
that'consists bf the'wake or approximately the remaining 75 to 80 percent of the boundary 
layer. In the inner region, a mixing length formulation for the eddy viscosity is used 


dinner == 




( 2 ) 


where 


£ = ky 


1 



(3) 
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with 


A=aJ-±£= 

k = 0,4, and = 26. In the outer region the eddy viscosity is represented by 

Wr = 0-0168 ^ 


where 



and 


1=1 + 




(4) 


(5) 


( 6 ) 


(7) 


In the calculations, (inner used in the boundary layer until y is large enough that 
^iniier ^ ^uter> which point and beyond (outer used. The heat flux vector is 
approximated by 


q 



p( 

Pr™ 


CpVT 


( 8 ) 


with Pr = 0.72 and Prrp = 0.9. The model described in equations (2) to (8) is denoted 
in the remainder of this paper as the "baseline model." 

In addition to the baseline model, various other models or variations of the base- 
line model were investigated. Calculations using these various models will be described 
in the "Results and Discussion" section. 

The calculations were started at a streamwise station (x - X[)/5o = -5.47 using 
the experimental turbulent boundary-layer velocity profiles and other information, such 
as shock -wave location, that is consistent with the experiment. To indicate the resolu- 
tion of the near -wall region of the boundary layer achieved in the numerical solution, a 
law of the wall velocity profile is shown in figure 2. The symbols are at the computa- 
tional mesh points and, as can be seen, at least one point is in the laminar sublayer 
where u+ = y+. The flow -field quantities at the mesh points were found by interpolating 
experimental data. Below the region where experimental data were available (y"*" less 
than about 500), the quantities were fo\ind by setting r = t^, neglecting streamwise 
derivatives in the numerical solution using the baseline turbulence model, and constrain- 
ing the velocity profile to pass through the experimental data beyond y"*" = 500. The skin 
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friction coefficient Cj = 1.30 x 10“^ calculated in this manner is in close agreement 
with the experimentally determined value obtained using a Preston tube (ref. 10) and 
with that obtained from examination of the velocity profile (ref. 12). 

RESULTS AND DISCUSSION 


ft - Comparison of Baseline Model With E?q)erimental Results . , 

The e:q)erimental results of reference 10 are compared with calculations using the 
baseline model in figures 3 to 6. Figure 3 contains a comparison of. experimental and, 
calculated wall pressures as a fimction. of the streamwise coordirute employed in ref- 
ereii'ce 10, namely - (x -Xi)/6o. Agreement between experiment and theory using the 
baseline model is poor.' Tn figure 4 the computed skin -friction coefficient is compared, 
with that obtained from velocity profiles from .reference .10 analyzed by. the method of . , 
reference 12. ^ The separation and reattachment points deduced- from. surface oil -flow ' 
observations (ref. 10) are also indicated. Not only is the. forward extent of the separa-. 
tio'n underestimated by the baseline -model calculation but the. rearward extent is overrc 
estimated. As a result of the' apparent inability of the baseline model to describe the , 
observed variation of pressure and skin-friction coefficient, various modifications to 
that model were made. 


Pressure. Gradient, Modification . 


To allow for effects of pressure gradients, Cebeci (ref. 13) recommends (in effect) 
replacement of A^ = 26 in equation (4) with 


A^ = 26 


+ 11.8 


15 

Me dp/dxV 
^e ^e 


( 9 ) 


where in the present study t = + y(ap/ax) (ref. 13) while the eddy viscosity in the 

outer flow remains 'unchanged^ According to Cebeci, the above correction ’factor is ' 
applicable in the presence of favorable as well as imfavorable pressure gradients. In 
the present application the modification increases the eddy viscosity in regions ahead;- 
of the separation as well as aft of reattachment. 


The effect of the Cebeci modification on skin friction is shown in figure 5. The 
separation bubble is moved slightly forward in the direction of the experiment, but the 
bubble is also extended farther rearward worsening the agreement with the measure- 
ments. Clearly some modification other than that of Cebeci (or in addition to) is needed 
for the class of separated flows considered here. 
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Escudier Formulation 


An alternative to equation, (5) of the baseline model that has been used for attached 
flows (ref. 14) is 

f 4 . = c2 3u . 9v 

‘^outer ^max ay ■ ax 

with fmax ~ 0-095. Results of using the Escudier model are showii also in figure 5. 

The naodel, however, is difficult to ^ply in the present calculation since the local width 
of the boundary layer 6 cannot be defined accurately enough in the numerical calcula- 
tion because of velocity gradients introduced by curvature of .the incident, induced,, and 
reflected shock waves. Slight clmnges in the definition of 6 . can change its value by .. ^ 
more than a-factor of 2. However, as observed in the experiment,: the values of j 6. do. 
not change significantly from the upstream ;value ahead of the separation. Accordingly, 
jfmax on the iqpstream boundary-layer thickness ^f^iax ~ 0.095oj has been used in 

several calculations to be described. Figure 6 shows calculated profiles of. pe in the 
separated region according to the baseline model and the Escudier formula. Although, 
as is well known, the two formulations produce similar variations for attached flows with 
small pressure gradients, they appear to produce greatly different variations for the 
present separated flow. 

Applying the Escudier formulation in the present calculations greatly reduced the 
extent of the separation along the wall. As evident in figure 5, the prediction of upstream 
influence was substantially worse than that of the baseline model, although the rearward 
extent of the bubble agreed better with indications from oil flow in the experiment. 

Despite the poor overall agreement shown, the Escudier model can be used as a guide in 
modifying the turbulence model. The better agreement with reattachment location indi- 
cates one type of modification to the turbulence model needed to obtain better agreement 
with the experimental data. Namely, relative to the baseline model, larger values of 
eddy viscosity are needed near reattachment to confine the rearward extent of the bubble. 
However, smaller values are probably needed near separation to promote the proper 
\q)stream influence. A model is now discussed that, in fact, does possess these 
characteristics. - 


Modification to Account for Noneqiiilibrium Effects 

For separated supersonic flow over a compression corner Shang and Hankey (ref. 8) 
have obtained remarkably good agreement with experiment by use of a simple modifica- 
tion of an algebraic turbulence model of the form (in the present notation) 


pe = (pe)o + (P^)equilibrium ■ (P^)o 


1 /1 _ AYn 

-a(x - xo) 1 

3 Y 

L 6o jj 


(10) 
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with a = 1/10. This type of modification (with a = 1/10) has been proposed by many 
authors arid has been shown by Rose and Johnso^n (ref. 11) to agree generally with the 
values of £max deduced from direct measurements of the Reynolds shear stress. A 
value of 0 ! between 1/5 and 1/15 is consistent with e}q)erimental observations, as noted 
in reference 11. 

AppUcation of equatira (10) in the present problem with a = 1/10 led to resiilts 
(not shown) simil^ to those obtained from the Cebeci viscous sublayer jmodification. 

Namely (relative to the baseline model), the upstream influence was inaproved, but the 
bubble was also extended far downstream in gross dis^reement with the experiment.. 

The observed effects of various turbulence models on the calculations (^scribed , n 
so far led to the following attempt? to modeb the turbulent shear stress. A formula like ; - 
equation '(10) should apply only to" the flow, outside the- separation bubble since there, probr iw 
ably shovild be no upstream history passing through the front of the bubble (i.e., the tur-. 
bulence within the bubble is in local equilibrium). In the present study, equilibrium 
within the bubble was insured by applyirig equatiori. (10) only if (pc) equilibrium exceeded 
(P€)q which does, in fact, allow a local equilibrium condition to prevail in the separation 
bubble. Since the Escudier formula produces the large values of eddy viscosity needed 
to limit the rearward extent of the bubble, (p€)equQii 3 j.imn is determined from that 
formulation. For the same reason in the present case, a value of a = 1/5 . was found to 
be better than a = 1/10. Iri other words the baseline model is modified by the use of , 
equation (10) with 

if (P^) equilibrium (P^)o 

a = < (11) 

(1/5 otherwise 

where (pe) equilibrium is determined from 


(P^outer) equilibrium “ 


^max 


3u ■ 9V 
9y 9X 


^max - 0.096 q 

(^^inner) equilibrium equations (2) to (4) of the baseline model if £ < fmax* 

Profiles of pe produced by these modifications are compared with those from the base- 
line model in figure 7. The desired variation in the values of pe is obtained; that is, 
upstream near separation pe is lower than the baseline model and downstream near 
reattachment pe is higher. 

In figure 8, calculated wall pressure distributions are compared to the experimental 
measurements. The upstream influence predicted by the relaxation model is in better 
agreement with the measurements than the prediction from the baseline or Escudier 
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models. However, the wall shear stresses shown in figure 9 indicate that prediction of 
the reattachment point is still poor. Velocity profiles near the front of the bubble pre- 
dicted by the baseline, relaxation, and Escudier models are compared with measure- 
ments in figure 10. Again, the results using the relaxation model are in much better 
agreement with the data than the simpler models. 

Finally, locations of the experimental shock waves and expansion from a tracing 
of a schlie'ren photograph of "tte flow field are compared in figure 11 with the locations 
of the shocks and expansion predicted using the relaxation model. The good agree- 
ment shown for the relaxation model would' not be obtained by any other turbulence 
models investigated in the present study since none of the other models predict the-. 
proper'^upstream influence.' One can thus conclude that the flow-field properties, such, 
as re^ons'of shock waves and expansions, can be adequately predicted when a; valid .• 
turbulence model is used. . ; . 

CONCLUDING' REMARKS ' , 

■ The complete, time -averaged conservation equations have been numerically solved 
using several algebraic turbulence models. Information has been gained on the effects 
of the various turbulence models applied to a separating and reattaching compressible 
turbulent boimdary' layer. The information has led to the formulation of a relaxation 
model for the Reynolds shear stress that produces reasonable agreement with a set of , - 
experimental data. 
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Figure 1.- Schematic drawing of flow field to be calculated. 


O UPSTREAM VELOCITY PROFILE 
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Figure 2.- Upstream velocity profile in law of the wall coordinates 
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Figure 3.- Comparison of measured pressure with calculation. 



Figure 4.- Comparison of predicted and experimental skln-frlctlon coefficients 

using baseline model. 
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O EXPERIMENTAL VALUES PROM 
VELOCITY PROFILES M. 13) 

■ SEPARATION AND REATTACNMENT 
: POINTS FROM OIL FLOW Iraf. 10) 

CALCULATIONS 



Figure 5.- Comparison of predicted and experimental akln-frictlon coefficients 

using various models. 

x-Xj/6^ = 0.097 



Figure 6.-*. Comparison of Escudler and baseline eddy-vlscoslty profiles 

In separation region. 
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Figure 9.- Comparison of predicted and experimental wall shear stress 

using various models. 


X-Xi 


= - 1.597 



U, m/sec 


Figure 10.- Comparison of predicted and experimental velocity profiles 

in separation region. 
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PREDICTED EXPERIMENTAL 

O INCIDENT SHOCK ESSJ (FROM SCHtlEREN 

0 INDUCED SHOCK' E2S] PHOTOGRAPH OF 

Cl EXPANSION ^ REF. 10) 

A REFLECTED SHOCK 



Figure 11,- Comparison of predicted and experimental flow fields 

using 'relaxation model. 



